Lossless Negative Refraction in An Active Dense Gas of Atoms 
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Current designs for negative index materials are passive, such that they suffer from significant 
losses that prohibit applications, in particular towards optical frequencies. We predict, based on 
nonlinear optical Bloch equations, negative refraction with zero absorption in an active, i.e., am- 
plifying, dense gas of atoms. External laser fields are used to tune the medium's electromagnetic 
response and an additional weak pumping field controls its absorption properties. Metastable Neon 
is identified as a suitable experimental candidate at infrared frequencies to implement a lossless 
active negative index material. 



Over the past years, tremendous progress has been 
accomplished in the field of negative refraction mostly 
related to the ongoing miniaturization of artificial struc- 
tures [l|,l2l- In such metamaterials, negative refraction is 
accompanied by a substantial amount of absorption espe- 
cially towards higher frequencies. These losses typically 
occur since the refractive index becomes negative only 
close to electromagnetic resonances where the absorption 
is high. All known metamaterials are passive devices and 
recent proposals to include gain materials to compensate 
for losses have not yet been successfully realized [3]]. As 
a result, the relevant figure of merit, the ratio between 
the real and the imaginary part of the refractive index 
FOM — |Re(n)/Im(n)| for high-frequency metamaterials 
is currently only of the order unityj)rohibiting most of 
the possible applications [J, [sj, [g, [Zl, |8|] . As a different ap- 
proach, recently, negative refraction without metamateri- 
als has been proposed in dense gases of atoms [9|, llO|, [ll | . 
Gases naturally have a macroscopic extend in all spatial 
dimensions, unlike high-frequency metamaterials which 
typically are produced layer by layer on a surface [12| . 
The proposed systems, however, have in common with 
current metamaterials that they are also passive and 
lossy. 

In this Letter, we demonstrate how truly lossless nega- 
tive refraction in an active, i.e., amplifying negative index 
medium can be realized. A dense gas of atoms is exposed 
to control laser fields that create coherence between dif- 
ferent internal atomic states. Quantum interference ef- 
fects reduce the absorption in the gas and an additional 
weak incoherent pumping field is used to render the sys- 
tem completely lossless: Im(n) = 0; or transfer it into 
an active, amplifying state where Im(n) < 0. Changing 
the intensity of the pumping field continously allows for 
a controlled transition from a passive to an active state. 

We employ nonlinear Bloch equations [13j| to account 
for dipole-dipole interactions, which goes beyond stan- 



dard local field corrections used in this context so far [11 1 . 
The feasibility of our scheme is demonstrated for the case 
of metastable Neon, where we find negative refraction in 
the infrared range at a wavelength of about A — 5/im. 
The required energy level scheme (see Fig.H]), could also 
be realized in other atomic species or solid state systems 



like doped semiconductors or quantum dot arrays. 

Negative refraction occurs if both the electric and the 
magnetic component of an electromagnetic probe wave 
couple near-resonantly to the system. Thus, in atomic 
gases, near-degenerate electric and magnetic dipole tran- 
sitions are required. The level structure of metastable 
Neon features such a pair of near-degenerate dipole tran- 
sitions, see Fig. [T] [10]. Additional laser fields with elec- 
tric component E^ and Ec greatly enhance the magnetic 
response, since they prepare a coherence on the mag- 
netic probe transition. We are interested in the linear 
response of the medium to the weak probe field of fre- 
quency Ljf) with electric component Ef, and magnetic field 
component H(,. The induced electric polarization P and 
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FIG. 1: (Color online) Five-level scheme with probe field E;,, 
Bt, (red) and coupling fields Ea, Ec, E^ (blue) of frequen- 
cies uji,, uia, iOc and ujd, respectively. The electronic states in 
jL-coupling notation refer to the special case of metastable 
Neon [iJ], with corresponding wavelengths A;, — 5.4^m, 
\a — 704nm, A^ — 352nm and A^ — 1.05/im. The magnetic 
and electric probe field transitions |1) — |2) (Ml) and |3) — 14) 
(El) are energetically degenerate up to a small energy gap 
A = 3.33 cm~^. An incoherent light field (black) acts as a 
pump rate r34, r43 between states |3) and i4). 



magnetization M at the frequency ujh are given by [H 

P = XEE{'^b)'^b + ^Bff (w6)Hb/47r 
M = XffHK)Hb + CHBK)Ef,/47r , 



(1) 



where XeEi Xhh are the electric and magnetic suscepti- 
bihtes, respectively, and S,EH^ Che are so-called chirality 
coefficients that describe the cross-coupling between elec- 
tric and magnetic fields in a chiral medium. The refrac- 
tive index depends on the probe wave polarization and 
we obtain for circular cr -polarization 
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where we have defined the electric permittivity e = 1 + 
'^t^Xee and the magnetic permeability /i = 1 -I- ^-Tixhh 
in the usual way. 

To find the linear response coefficients Xaa and ^q/j 
(a ^ /3 G {i?, iJ}), we consider the time-evolution of 
a single five-level atom with density matrix p which is 
governed by the Master equation [16|] 



^ = ^[^"^]-E^(b')01'°+'°b')01-2|fc)01Pb->(^i), 

(3) 



i,fc 



with the system Hamiltonian given by 
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(4) 



Here, Cj is the energy of state |j) and the other terms 
describe the interaction with the laser fields in the long- 
wavelength and dipole approximations. The fields have 
frequencies Wi with i G {a, 6, c, d}. The laser detuning 
on transition \i) — |j) reads simply M^j — hujij + eij, 



and eij accordingly, if we intro- 



where uJij = uji ^ 

duce the following combination of angular frequencies 

LUi — LOd, L05 — 0. The electric control field compo- 



nents E, 
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c.c. with polarization vector e^ 



{i e {a,c,d}) give rise to the complex Rabi frequencies 

rial = dgi • Ea/ft, ^42 = d42 • Ec/K ^^54 = d54 • Ed/H, 

where d^ — (z|d|j) is the electric dipole operator be- 
tween states \i) and |j). Note that the Hamiltonian 
Eq. ^ contains instead of the external probe fields Eb, 
H(, the actual (local) fields inside the medium E^ = 



Eft 



^P, H 



L — 



Hft 



^M 17]. This accounts for 
dipole-dipole interactions which cannot be neglected for 
large atomic densities. 

The second part of the master equation describes the 
spontaneous decay of excited states and the decoherence 
that arises due to elastic collisions. The decay rate on 



transition \i) — \j) is denoted by jij and set to 7^ = 7, if 
|z) — |j) is electric dipole allowed (El) and set to -fij = a^7 
for metastable E2, Ml transitions, where a = 1/137 
is the fine- structure constant. A typical decay rate in 
metastable Neon in angular frequencies is 7 = 10^—. To 
include the effect of elastic collisions, we add a coUisional 
decay constant jc to the decay rates of the off-diagonal 
density matrix elements: ^xy = Y^jilxj + lyj)/'^ + Ic 
(x ^ y e {1, . . . , 5}). In the following, we set 7c = 7- 
A numerical analysis shows that our system is quite ro- 
bust against changes in the decay rates. In particular, our 
model does not rely on the preservation of an electromag- 
netically induced transparency ground state coherence as 
in previous models [111 ] , which is challenging in particular 
at high particle densities. 

The steady-state coherences P43 and p2i determine 
the induced polarization and magnetization at the probe 
field frequency as P = iV (p43d34 -I- c.c.) and M = 
N (p2iAti2 + c.c), where N is the density of atoms, and 
/^i2 — (1|m|2) is the magnetic dipole moment of the 
Ml-transition. It turns out that large atom densities 
are required to obtain negative refraction, meaning that 
many particles are found in a cubic wavelength volume 
N\^ ;2> 1. Under such conditions, the local probe fields 
El, Hi inside the medium may considerably deviate 
from the externally applied probe fields Eft, Hft in free 
space, since they contain contributions from neighboring 
atoms. The fields are related by the Lorcntz-Lorenz (LL) 
formulas [11]: E^ = Eft -H f^P and Hl = Hft + ^M. 
Different approaches have been proposed in the literature 
to solve this problem. We found the simplest approach, 
which is expected to be valid for moderate densities and 
consists of expanding the steady-state density matrix el- 
ements P21 and P43 to linear order in the local fields E/,, 
Hi and only then substitute the external fields for the 
local fields using the LL- relations, to be insufficient to 
describe our system. Therefore, we relate the local fields 
via the LL-relations to the external fields already in the 
density matrix equations of motion 13| and obtain a set 
of non- linear equations. This has the advantage that in 
the framework of the LL-formulas the local field effects 
are treated without approximations. Since the equations 
are non-linear, however, one typically requires a numeri- 
cal analysis. 

We numerically integrate the nonlinear equations un- 
til the system has reached its steady state for (twenty) 
different electric field amplitudes Eb, holding the mag- 
netic field amplitude H}, fixed. Linear regression of the 
relevant coherences p2i , P43 as a function of Ei, allows to 
extract the response coefficients as slope and y-axis inter- 
cept from Eqs. ([T|). We note that in the case of a circu- 
larly polarized probe field, these equations describe scalar 
relations, since H^ — B;, = =F*Eft for (T*-polarization. In 
particular for the electric polarization 
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where P43 = p43e"^'' is the slowly varying part of the 
coherence, we obtain for a cr" -polarized probe field 



P = 2iVp43(i34 = xee 
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with the small expansion parameters we — di^Et/h'y and 
wb = l^2iBb/fi'y- The dipole moments dij and /i^ can 
be evaluated via the respective spontaneous decay rates 

as J3jijhc^/{Aujf,). By hnear regression of P{we) — 

mpWE + bp one finds 

^43 mp 
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Analogously, one obtains the response coefficients for the 
magnetization M = tumwe + bi^i: 
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It turns out that the chirality coefficients depend on the 
relative phase of the applied laser fields. Since they are 
typically several orders of magnitude smaller than the 
direct coefficients e and /i, however, they do not con- 
tribute significantly to the index of refraction. Still, in 
our calculations we average over this phase to simulate 
an experiment without phase control, which effectively 
reduces the magnitude of ^eh and S,he almost to zero. 

In our setup, the coupling laser fields serve different 
purposes: the fields E^ and Ec enhance the magnetic 
response by inducing the coherence p2i which drives the 
magnetic dipole moment of the atom. We choose the 
frequencies of both fields to be equal uJa — ^c and almost 
resonant with the transition 1 1} — 1 3): (Jsi ^ fisi [18[. The 
third laser field E^ is resonant with the transition |4) — 
1 5) (1^54 = 0), causing an Autler-Townes-splitting of the 
fourth level into two dressed states at £4 — > ei±h\^^4\/2. 
Thereby, electric and magnetic probe field resonances are 
shifted with respect to each other. For zero field, they 
are separated in Neon by the gap A = (£43 — e2i)/h — 
3.33 cm-i = lOOGHz = 27r 10^*7, since 7 = 10^/(27r) Hz 
is a typical value for the spontaneous decay rate in Neon. 
Note that in angular frequencies this corresponds to 7 = 
2Q7r|d jj^ order to close the gap, one has to apply a 
strong laser field with ^^4 « 10^7. Yet we find that 
neighboring transitions are still detuned at the least by 
15ri54, making such large light shifts feasible without 
inducing unwanted transitions. Thus the field O54 allows 
to circumvent the problem of non-degenerate electric and 
magnetic probe field transition 9\. It is important to 
note, however, that it is not mandatory to close this gap 
completely, because the states |2) and |3) are connected 
by a two-photon transition induced by the fields Ec and 
Eb which becomes important for nonzero gap frequencies 
A. 

This two-photon transition is the motivation to apply 
in addition an incoherent light field resonant with the 
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FIG. 2: (Color online) Electric permittivity e and magnetic 
permeability /i as a function of 521 for different pump rates 
n = 0, ra = 0.2512, rg ^ 1, n = 1.679, rg = 1.698, 
rg — 1.799 (in units of 10~^7) between states |3) and |4). 
The arrow denotes the direction of increased pumping. Other 
parameters are set to ^42 = 5.6, Q31 — 6.3 • 10~^, A' = 560, 
S31 — — 1Q~^, S54 = (aU in units of 7). The probe field 
wavelength is A = 5/im and the density A'^ — 2.5 ■ lO'^^cm""^. 



transition |3) — |4), which transfers population between 
the levels with rates r34 = r43 = r changing the de- 
cay rates from 743 -^ 743 -I- r43 and accordingly for 734. 
Together with spontaneous emission from |4) to |2), it 
effectively pumps population into |2). One finds that for 
P22 > P33 the probe field is amplified by means of the 
two-photon process, since it is more likely for the tran- 
sition to occur in the direction |2) -^ |3). This direction 
involves the emission of a probe field photon and we ob- 
serve gain in the electric probe field component. Since 
we choose the laser frequencies Ua — i-^c, this two-photon 
resonance is always located at the position of the mag- 
netic probe field resonance which occurs at S21 = 0. The 
two-photon virtual internrediate level is separated from 
|4) by the effective gap A' = —A -|- $754/2. If we choose 
A' sufficiently large, the two-photon transition will be 
the dominant process around zero probe detuning. 

In Figure [21 the permittivity and the permeability of 
the medium are shown as a function of the probe field 
detuning from the magnetic resonance 621. The resulting 
refractive index can be seen in Fig. [3l Without incoher- 
ent pumping the system is passive and thus e, /i and n 
have a positive imaginary part, i.e., the medium absorbs. 
Whereas absorption is small for the electric response due 
to local field effects, the losses are significant in the mag- 
netic component. We observe power broadening in all 
response functions due to the nonlinearities that appear 
in the master equation to include near-dipole-dipole ef- 
fects. Gradually increasing the incoherent pumping rate, 
we find that the imaginary part of e around zero detun- 
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FIG. 3: (Color online) Refractive index n for different pump 
rates as a function of ^21 • The parameters are the same as in 
Fig. [2I Arrows indicate the direction of increased pumping. 
Without pumping, the imaginary part is strictly positive, as 
expected for a passive system. On increasing the pump rate, 
the system turns active, and negative refraction is obtained 
without absorption or even with amplification. 
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sive and shows significant absorption (FOM < 1). In- 
creasing the pump rate, the imaginary part of e becomes 
negative which compensates for the losses via the mag- 
netic component. The imaginary part of the refractive 
index decreases (FOM ^ 1) and finally vanishes, while 
the real part remains negative. Thus negative refraction 
with a large figure of merit can even occur for positive 
Re(e, /i) if |Im(e, ;u)| ^ and the losses via the magnetic 
field component are compensated by gain via the electric 
component. Since we can continously change the pump 
rate, we can identify the physically correct branch if we 
follow the path in the complex plane starting from the 
root with positive imaginary part for the passive system. 

To summarize, we have predicted negative refraction 
with truly zero loss in a dense gas of atoms. Employing 
the transition to an active medium by means of an inco- 
herent pumping field allows to study a qualitatively new 
parameter range not accessible with current devices. Due 
to its wide tunability and control, our system is not only 
interesting from a proof-of-principle point of view, but 
also promises the enhancement of optical effects which 
are severely degraded by losses accompanying negative 
refraction in current metamaterials. 



FIG. 4: (Color online) (a) Lossless negative refraction for 
a pump rate of r = 1.718 • 10~^ 7. Re(n) (solid red line), 
Im(n) (dashed blue line), figure of merit FOM=|Re(n)/Im(n)| 



(solid black line), (b) Paths of e, fi, n , and ±\/n^ in the 
complex plane parameterized by the pumping rate r. Step 1 
corresponds to zero pumping ri, the system is in a passive 
state and the physical solution is denoted by a filled circle 
(red). Increasing the pump rate continuously, we follow the 
path indicated by the arrow. Step 2 is at an intermediate 
pump rate in Fig. [S] where the system is still absorptive. 
Finally, step 3 corresponds to lossless negative refraction. 



ing turns negative, indicating that the two-photon tran- 
sition |2) — |3) amplifies the probe field. The magnetic 
permeability is mostly unaffected by the incoherent field. 
The real part of n is negative for all detunings, which 
even includes regions of positive Re(/i), but absorption is 
dominant there. Without pumping, the imaginary part 
is strictly positive, as expected for a passive system. 

Increasing the pumprate r decreases the absorption 
over the whole range of displayed probe detunings. Fi- 
nally, negative refraction occurs without absorbtion or 
even with amplification, see Fig. |3](a). In this frequency 
range, the FOM becomes very large due to the vanish- 
ing imaginary part. The crossover from passive to ac- 
tive induced by the incoherent pumping can easily be 
understood by following the path of e, /i, 71^ and n in 
the complex plane as a function of the pump rate r, see 
Fig. mb). Without pumping (step 1), the system is pas- 



[1] V. G. Veselago and E. E. Narimanov, Nature Materials 

5, 759 (2006). 
[2] V. M. Shalaev, Nature Photonics 1, 41 (2007). 
[3] A. D. Boardman, Y. G. Rapoport, N. King, and V. N. 

Malnev, J. Opt. Soc. Am. B 24, A53 (2007). 
[4] R. A. Shelby, D. R. Smith, and S. Schuhz, Science 292, 

77 (2001). 
[5] H. J. Lezec, J. A. Dionne, and H. A. Atwater, Science 

316, 430 (2007). 
[6] E. E. Narimanov and V. A. Podolsky, Opt. Lett. 30, 75 

(2005). 
[7] D. R. Smith, D. Schurig, M. Rosenbluth, S. Schuhz, S. A. 
Ramakrishna, and J. B. Pendry, Appl. Phys. Lett. 82, 
1506 (2003). 
[8] R. Merlin, Appl. Phys. Lett. 84, 1290 (2004). 
[9] M. O. Oktel and O. E. Miistecaplioglu, Phys. Rev. A 70, 
053806 (2004). 
[10] Q. Thommen and P. Mandel, Phys. Rev. Lett. 96, 053601 

(2006). 
[11] J. Kastel, M. Fleischhauer, S. F. Yelin, and R. L. 

Walsworth, Phys. Rev. Lett. 99, 073602 (2007). 
[12] G. Dolling, M. Wegener, and S. Linden, Opt. Lett. 32, 

551 (2007). 
[13] C. M. Bowden and J. P. Bowling, Phys. Rev. A 47, 1247 

(1993). 
[14] http://physics.nist.gov/Pubs/AtSpec/; 
http://physics.nist.gov/PhysRefData7] 
[15] J. B. Pendry, Science 306, 1353 (200l). 
[16] M. O. Scully and M. S. Zubairy, Quantum Optics (Cam- 
bridge, 1997). 
[17] J. D. Jackson, Classical Electrodynamics (John Wiley & 

Sons, 1998). 
[18] M. Mahmoudi and J. Evers, Phys. Rev. A 74, 063827 



(2006). 



